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Abstract 

The reorientation transition of the magnetization of ferromagnetic 
films is studied on a microscopic basis within a Heisenberg spin model. 
Using a modified mean field formulation it is possible to calculate 
properties of magnetic thin films with non-integer thicknesses. This 
is especially important for the reorientation transition in Ni/ Cu(OOl), 
as there the magnetic properties are a sensitive function of the film 
thickness. Detailed phase diagrams in the thickness-temperature plane 
are calculated using experimental parameters and are compared with 
experimental measurements by Baberschke and Farle (J. Appl. Phys. 
81, 5038 (1997)). 

The direction of the magnetization of thin ferromagnetic films depends 
on various anisotropic energy contributions like surface anisotropy fields, 
dipole interaction, and eventually anisotropy fields in the inner layers. 
These competing effects may lead to a film thickness and temperature 
driven reorientation transition (RT) from an out-of-plane ordered state 
at low temperatures to an in-plane ordered state at high temperatures at 
appropriate chosen film thicknesses. Experimentally, this transition has 
been studied in detail for various ultra-thin magnetic films [[Tp. Recently, 
it was found by Farle et al. [0] that ultra-thin Ni-films grown on Cu(001) 
show the inverse behavior: the magnetization is oriented in-plane for thin 
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films and at low temperatures and perpendicular for thick films and at 
high temperatures. 

It has been shown by various authors [^)-§P that the mechanism respon- 
sible for the temperature driven transition can be understood within the 
framework of statistical spin models. The RT from an out-of-plane state 
at low temperatures to an in-plane state at high temperature is found to 
be due to a competition of a positive surface anisotropy and the dipole 
interaction. It can occur because the surface anisotropy has a different 
temperature dependence that the dipole exchange anisotropy and van- 
ishes more quickly when approaching the Curie temperature The 
thickness driven reversed RT in ultra-thin Ni-films is argued to have its 
origin in a stress-induced positive uniaxial anisotropy energy in the inner 
layers with its easy axis perpendicular to the film [|2|]. This anisotropy 
is in competition with the dipole interaction and with a negative surface 
anisotropy. This scenario can indeed lead to a reversed temperature driven 
RT, but now the reduced surface magnetization plays a crucial role 
A third type of RT has recently been found theoretically, where the magne- 
tization switches from perpendicular to in-plane direction with increasing 
temperature, but with decreasing film thickness [§]. 

In this paper we will focus on the phase diagram of Ni/Cu(001) and 
on the influence of microscopic fourth-order anisotropy terms. The calcu- 
lations are done in the framework of a classical ferromagnetic Heisenberg 
model consisting of L two-dimensional layers on a face centered cubic 
(001) lattice. The Hamiltonian reads 

n = - 3 - E^-^-E^(^) 2 +^(^) 4 

(ij) i 

+ ^ ' *3 ~ 3r ij 5 & ' TiiWij ■ Sj), (1) 

/ ij 

where Sj = (sf, s\ , sf) are spin vectors of unit length at position fj = 
iff , rf , rf ), and = r { — fj. J is the nearest-neighbor exchange coupling 
constant. The uniaxial and fourth-order anisotropics D l n are position- 
dependent: D l n = if spin i is in the top layer, D l n = D% if spin i is 
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in the bottom layer, and D l n = otherwise. Finally, cu = fx fj, 2 / Ana 3 is 
the strength of the long range dipole interaction on a lattice with next- 
neighbor distance a (/x is the magnetic permeability and /i is the effective 
magnetic moment of one spin). 

The Hamiltonian Eq. (P is handled in a molecular-field approxima- 
tion [§]. In the following we assume translational invariance within the 
layers and furthermore that the magnetization is homogeneous inside the 
film and only deviates at the surfaces Therefore we can set (sj) = (s A ) 
if Si is a spin in a volume layer (A = V) or in the top or bottom surface 
layer (A = S, S'). We will focus on the case L > 2, i.e. both surfaces are 
complete. The resulting system with three mean field spins has the effec- 
tive interactions x XfM , where Xg is either the number of next neighbors, zg, or 
dipole sum in fcc(OOl) geometry, $ 5 = , respectively, 

between layers A and A +5. We have xss' = #s's = 0, xss = xs'S' = %q, ^sv = 
^S'v = %i/ ^vv = + 2xi(l - Ly 1 ), and x ys = ^vs' = XiLy 1 , with L s = 
L$> = 1 and L v = L — 2. The constants zg and Qg are z = z\ = 4, <3> — 9.034 
and $i = 1.429. ^s>i can be neglected in our approach. Note that the cou- 
pling between the volume spin and the surface spins is asymmetric. With 
this method we are not restricted to films with integer values of the thick- 
ness anymore. The effective interactions x\^ enter the mean field Hamil- 
tonian via the mean fields in layer A, h\ = X^ =1 Jz\^m^ + I^Wm^, with 
the dipole exchange anisotropy W = diag(l/2, 1/2,-1). In order to deter- 
mine the phase diagram of this Hamiltonian in the thickness-temperature 
plane, we directly calculate the phase boundaries using a stability analy- 
sis of the mean field free energy. Defining the function e m ^ n = evmin(A) 
which returns the smallest eigenvalue e m [ n of the matrix A, we can cal- 
culate the effective anisotropy K 2 from the Hessian of JF MF with respect to 
the azimuth angle i3x of the magnetization m A [ ]T0| ] 

)■ (2) 



Now consider the phase with magnetization parallel to the film normal 
z, where -Q\ = 0. At thicknesses L and temperatures T where this phase 
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is stable, we have K 2 (L, T, 0) > 0. At the phase boundaries both to the 
canted phase and to the paramagnetic phase, K 2 (L, T, 0) becomes zero, 
while it is negative in the canted phase and in the phase with in-plane 
magnetization. On the other hand, K 2 (L,T,ir/2) is positive when the 
in-plane phase is stable and becomes negative in the phases with per- 
pendicular component of the magnetization. In the paramagnetic phase 
K 2 (L, T, = for all Hence the reorientation temperatures are given 
by K 2 (L, T? y (L), 0) = and K 2 (L, T r z (L), § ) = 0. Next we will examine the 
phase diagram of this model for parameters measured on Ni/Cu(001) 0|: 
The exchange interaction is approximated via the Curie temperature T c = 
631 K of bulk Ni, using the mean field formula 3T C = z J with Zf cc = 12 
for the isotropic classical Heisenberg model to give J ~ 13.6 meV. This 
is a rather rough estimate, but note that the exact value of J has nearly 
no influence on the following results, as long as J is large compared to 

o 

the other energies in the model. With /i = 0.62/iB and a = 2. 5 A we 
get uj = 1.3 fieV ~ 10~ 4 J for the dipole constant. The uniaxial surface 
anisotropy is (-Df + Df)/2 = —60 fieV ~ —47a; and the uniaxial volume 
anisotropy is D% = +40 /zeV ~ 3 lw. Note that our surface anisotropics 
are the sum of the experimental surface and volume part, as our surface 
layers also count to the volume. For simplicity we assume that the bottom 
surface carries the volume anisotropics, hence we get Z)f ~ — 125u;. 

To check the influence of D% ±/ we show the phase diagram for the 
abovementioned parameter set with D% ± = 0. The phase boundaries 
are the thin lines in figure |1|. We obtain a reversed RT as expected for 
Ni/Cu(001), the reorientation thickness ranges from 5ML to 7ML which 
is slightly lower as in the experiment [p. We find that the canted phase 
exist, but it is rather narrow at these parameters. This canted phase is 
stable because we allow non-collinear magnetizations in different layers. 
To reproduce the experimental finding of the rather broad canted phase, 
we set the fourth-order anisotropics D% ± ~ —\D 2 The phase dia- 
gram of our model with these values and slightly modified D\ is depicted 
in figure [j] together with experimental data obtained by Baberschke and 
Farle [§]]. As the mean field theory always overestimates T c , especially for 
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low dimensional systems, the critical line T C (L) from [|p cannot be repro- 
duced well within our theory. Thus we normalized the temperature axis 
not with T c bulk , but with the Curie temperature of a film with L = 7 ML. 

In this paper we defined a modified mean field model to describe the 
RT in Ni/Cu(001). Using appropriate model parameters close to exper- 
imental findings, we find a rather nice agreement of theory and experi- 
ment. To reproduce the width of the phase with canted magnetization, we 
need to introduce microscopic fourth-order anisotropics D^±. 

This work was supported by the Deutsche Forschungsgemeinschaft 
through Sonderforschungsbereich 166. Discussion with K. Baberschke is 
gratefully acknowledged. 
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Figure 1: Calculated phase diagram for Ni/Cu(001). For the thick lines, 
the model parameters are u = 10~ 4 J, Df = -100a;, = Df = 24a;, 
Df ± = 24a;, D^ ± = Df ± = —4a;. The thick solid line represents the Curie 
temperature T C (L) of the film, the dotted and dashed lines are T r z (L) and 
T* y (L), respectively. The symbols are experimental measurements of the 
reorientation transition in Ni/Cu(001) taken from [|]]. The temperature 
axis is normalized to the Curie temperature of a film with thickness L = 
7 ML. The parameters for the thin lines are described in the text. 



